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Abstract 

The lattice dynamics of Ti02 in the rutile crystal structure was studied by a combination of 
thermal diffuse scattering, inelastic x-ray scattering and density functional perturbation theory. 
We experimentally confirm the existence of a predicted anomalous soft transverse acoustic mode 
with energy minimum atq = (1/2 1/2 1/4). The phonon energy landscape of this particular branch 
is reported and compared to the calculation. The harmonic calculation underestimates the phonon 
energies but despite this the shape of both the energy landscape and the scattering intensities are 
well reproduced. We find a significant temperature dependence in energy of this transverse acoustic 
mode over an extended region in reciprocal space which point to a strong role of anharmonicity 
in line with a substantially anharmonic mode potential-energy surface. The reported low energy 
branch is quite different from the ferroelectric mode that softens at the Brillouin zone center and 
may help explain anomalous convergence behavior in calculating Ti 02 surface properties and is 
potentially relevant for real behavior in Ti 02 thin films. 
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I. INTRODUCTION 


The transition metal oxide Ti02 with its mixed ionic and covalent bonding is a challeng¬ 
ing system in condensed matter physics with many interesting applications. Its low snrface 
energy with a particnlar energy landscape is, for example, very attractive for snrface chem¬ 
istry and photo-catalysis. In the rntile crystal strnctnre Ti02 is an incipient ferroelectric 
and shows a large and strongly temperature-dependent dielectric constant. The consequent 
high refractive index is valuable for a diverse range of technological applications such as di¬ 
electric mirrors, antireflecting coatings and pigments. Similarly to ferroelectric perovskites, 
the dielectric constant increases with decreasing temperature. The incipient ferroelectric 
behavior and the high dielectric constant are directly linked to the lattice dynamics and 
they were shown to be a consequence of a low-frequency transverse optical phonon mode 
A. 2 u^. Experimentally the phonon dispersion relation could be measured using neutron 
spectroscopy^. Effects of temperature and pressure on the static dielectric constants, fre¬ 
quencies of Raman active modes and Griineisen parameters were addressed by capacitance 
measurements and Raman scattering and a particularly strong dependence on both pressure 
and temperature was found for the A 2 U mode and related to structural and anharmonic 
contributions^’^. Ab initio lattice dynamics calculations using density functional perturba¬ 
tion theory (DEPT) were used to calculate phonon frequencies. Born effective charges and 
the dielectric permittivity^’'’^^. These calculations emphasize the importance of the mixed 
covalent ionic bonding of s and d orbitals which is linked to its large polarisability due to 
long range coulomb interactions between the ions. The large Born effective charges obtained 
could be explained by dynamical electron transfers during atomic displacements. Some of 
the phonon modes show a pronounced stress dependence: A significant softening of the fer¬ 
roelectric A 2 n phonon mode was found®. A recent theoretical work found an anomalously 
soft transverse acoustic (TA) mode quite separate from the A 2 U mode®. This mode shows 
a remarkable stress dependence in a region of reciprocal space unexplored by experiment 
and softens at around q = (1/2 1/2 1/4) under isotropic expansion. The atomic displace¬ 
ments of all atoms for this mode at this specific momentum transfer are perpendicular to 
the (110) surface. They allow for long ranged propagation of a displacement arising at the 
surface deep into the bulk. Under uniaxial strain the mode is predicted to soften for both 
expansion and compression. The phonon frequencies are very sensitive to slight changes in 
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the geometry and thus to the applied exchange-correlation functionals and accuracy of the 
pseudopotentials®. Very recently the importance of anharmonicity in the lattice dynamics 
of Ti 02 rutile has been addressed by inelastic neutron scattering and ab initio molecular dy¬ 
namics simulations®. Temperature dependent measurements of the phonon density of states 
revealed a signihcant increase in energy of the TA phonons with temperature. The phonon 
potentials of this mode were reported to develop remarkable quartic terms due to a dynamic 
change of the hybridization between Titanium and Oxygen atoms. 

Inspired by the work of Mitev et al.^ and motivated by the fact that soft modes are very 
sensitive to the approximations used in theory, we experimentally address the TA mode 
by inelastic and diffuse x-ray scattering and compare the results to harmonic DFPT cal¬ 
culations. We explore an extended region of reciprocal space by thermal diffuse scattering 
(TDS), measure the energy of the TA phonons with mode level resolution by temperature 
dependent inelastic x-ray scattering (IXS) and carefully compare the phonon energy land¬ 
scape and scattering intensities to the calculation. Finally we discuss details of the softening 
and importance of anharmonicity. 


II. METHODS 

A. Experiments 

Synthetic Ti02 single crystals from CrysTec GmbH (Berlin, Germany) grown by Verneuill 
process were used in this study. A rectangular bar with dimensions of approximately 0.2 x 
0.2 X 1.5 mm® was produced from a commercially available crystal by cutting and polishing. 
The crystal quality was checked by high-resolution x-ray diffraction and is visible in the 
reconstructed planes of diffuse scattering, see Section III. 

Single crystal IXS was carried out at beamline ID28 at the European Synchrotron Radia¬ 
tion Facility (ESRF) in Grenoble, France. The spectrometer was operated at 17.794 keV in¬ 
cident energy, providing an energy resolution of 3.0 meV full-width-half-maximum (FWHM). 
The momentum resolution was set to 0.3 (horizontal) x 0.9 (vertical) nm“^ FWHM. Energy 
transfer scans were performed in transmission geometry at constant momentum transfer (Q), 
selected by appropriate choices of scattering angle and sample orientation. Further details 
of the experimental setup can be found elsewhere^®. The temperature dependent study com- 
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bines data from experiments using a closed cycle Helium cryostat (20K, 50K, 75K, 122.5K, 
150K, 200K and 250K), a cryostream cooling system (105K) and no cooling device at all for 
the room temperature measurement. 

X-ray diffuse scattering was performed on beamline ID29 at the ESRF. Monochromatic 
X-rays with 16 keV energy were scattered form the single crystal in transmission geometry. 
The sample was incrementally rotated orthogonal to the beam direction by 0.1° over a 
full cycle while diffuse scattering patterns were collected in shutterless mode with a single 
photon counting pixel detector with no readout noise and pixel size of 172 pm (PILATUS 6M 
detector from Dectris AG., Baden, Switzerland). The orientation matrix was dehned using 
the CrysAlis software package (Agilent technologies, Santa Clara, USA). Reciprocal space 
reconstructions were prepared using locally developed software. Corrections for polarization 
of the incoming X-rays and for solid angle conversion associated with the planar projection 
were applied to the scattering intensities of individual pixels. 

B. Calculation 

Lattice dynamics calculations were performed in harmonic approximation using den¬ 
sity functional perturbation theory (DFPT)^^ as implemented in the CASTFP code^^’^^. 
In the calculations we employed local density approximation (LDA) with the exchange 
correlation functional by Perdew and Zunger^^. We used the plane-wave formalism and 
norm-conserving pseudopotentials of the optimised form^^. The LDA exchange correlation 
functional was chosen because the derived lattice parameters and vibrational properties are 
correctly described. Within the generalized-gradient approximations (PBF and PW91) the 
equilibrium crystal structure of Ti 02 is predicted to be unstable, which leads to imaginary 
phonon frequencies®’^®. We have used the pseudopotential from Mitev et al.^ which treats 
hve reference orbitals as valence states for Titanium and two for Oxygen. A plane wave 
cut-off of 1100 eV and an electronic grid sampling on a 4 x 4 x 6 Monkhorst-Pack grid en¬ 
sured convergence of forces to < 1 meV/A. The cell geometry was optimised employing the 
Broyden-Fletcher-Goldfarb-Shannon method^^ by varying lattice and internal parameters. 
For the cell parameters of the optimized cell we hnd a = b = 4.574 A and c = 2.944 A. 
Phonon frequencies and eigenvectors were computed on a 5 x 5 x 7 Monkhorst-Pack grid 
of the irreducible part of the Brillouin zone by perturbation calculations. Sum rules for the 
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acoustic branches as well as the charge neutrality at the level of Born effective charges were 
imposed with a maximum correction of 2.0 meV at F. The phonon dispersion throughout 
the hrst Brillouin zone was calculated using DFPT in conjunction with Fourier interpo¬ 
lation with a grid spacing of 0.005 in the cumulant scheme including all image force 
constants^®’^®. The maximum error in phonon energy from the interpolation is < 1.5 meV. 

Scattering intensities for inelastic and thermal diffuse scattering were calculated in hrst 
order approximation assuming the validity of both harmonic and adiabatic approximation. 
The formalism can be found elsewhere^°’^°’^^. Debye-Waller factors were calculated from a 
hne phonon grid with a spacing of 0.005 A“^. The momentum transfer dependence of the 
atomic scattering factors was taken into account using an analytic function with coefficients 
derived from Hartree-Fock wavefunctions^^. IXS spectra are resolved in phonon energy and 
the intensity is sensitive to the eigenvectors. TDS intensities are not resolved in energies, 
but due to their strong energy dependence, they are sensitive to low energy features. 


III. RESULTS 

Phonon dispersion relations along selected high symmetry directions in the Brillouin zone 
are shown in Figure 1. They combine the results from our calculation, IXS measurement 
at 295 K and previous INS data at 300 K^. INS and IXS data are in excellent agreement, 
as expected. Experimental phonon energies are reproduced by the calculation within 2.5 
meV or better, except for the transverse acoustic branch along T — Z and M — A, where 
the discrepancy rises to nearly 5 meV. Importantly, we experimentally conhrm the existence 
of the previously predicted low energy feature close to a reduced momentum transfer q = 
(1/2 1/2 1/4). IXS scans along M — A are summarized in an intensity map, see Figure 
2. The dispersion clearly shows dips at Q = (2.5 2.5 0.75) and (2.5 2.5 1.25) which both 
correspond to q = (1/2 1/2 1/4)^®. The shape of the dispersion is reasonably reproduced 
within the harmonic (OK) DFPT calculation in spite of a signihcant underestimation in 
phonon energies. 

The variation of the TA mode energy across the HHL plane is investigated by IXS and 
measured at a grid of Q-points on that plane. The extracted phonon energies of this partic¬ 
ular mode are presented on a color map and compared to the calculation in Figure 3. We 
note that the phonon energy landscape forms an elongated valley. The energy landscape is 
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FIG. 1. Phonon dispersion relation of Ti02 rutile, calculated by DFPT (full lines) and measured 
by IXS at 295 K (blue points) and INS at 300 K (red triangles)^. The position of q = (1/2 1/2 
1/4) is indicated by an arrow. 


nicely described by the calculation along [^^/]. It is also well reproduced along [hh^], except 
in the center of the image, where the valley appears to be a little broader (along [00.^]). A 
closer look at the absolute energies (note the different absolute color scales in the hgure) 
reveals that the calculated phonon energies differ by about 5 meV across the plane. Only 
the red areas have comparable phonon energies. 

The variation in 3D reciprocal space is investigated by TDS. Scattering intensities on 
reconstructed high symmetry reciprocal space sections and the 3D iso-intensity distribution 
around Q = (2.5 2.5 1) are shown in Figure 4 and compared to calculations. On the HKO 
plane we note that the intensity distribution of TDS, which contains contributions from 
all phonon modes, is well reproduced by the calculation across the plane. Indeed, in this 
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FIG. 2. IXS intensity of Ti02 crystal as a fnnction of energy transfer along Q = (2.5 2.5 0.5) to 
(2.5 2.5 1.5) (A — M — A), obtained from experiment at 295K (a) and calcnlation (b). Positions 
of q = (1/2 1/2 1/4) are indicated by arrows. The experimental map consists of eleven linearly 
spaced spectra with energy steps of 0.7 meV. The calculation was performed at 200 q-points. 

plane the calculated phonon dispersion is quite accurate. The intensity of diffuse scattering 
is also sensitive to the phonon eigenvectors. The excellent agreement gives high degree of 
conhdence that the calculation properly reproduces both energies and eigenvectors across 
this reciprocal space section. The diffuse scattering from the TA mode along T — M is 
particularly strong and forms streaks of high intensity for example between (0 2 0) and (2 0 
0). On the HHL plane one can identify the diffuse scattering from the soft TA mode. The 
shape of this feature is very similar to the phonon energy landscape (Figure 3). In fact, 
this mode is well separated from higher energy modes and contributes to 98 % to the total 
TDS at Q = (2.5 2.5 0.75). The effect of the calculation underestimating the phonon energy 
is clearly noticeable: Calculated diffuse scattering intensities are signihcantly stronger than 
experimentally observed. An illustration if the iso-intensity distribution gives insight to the 
energy landscape in 3D, see Figure 4 c and d. We note an elongated shape in HHL with 
very little variation orthogonal to it. 

The variation of selected phonon energies with temperature is shown in Figure 5. A 
signihcant temperature dependence of the soft TA mode at q = (1/2 1/2 1/4) and Z is visible. 
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FIG. 3. Phonon energy landscape of the soft TA mode on the HHL plane. The phonon energies 
depicted in color were extracted from IXS measurements taken at 105K on 90 Q-points (a) and 
calculated on 5150 Q-points (b). Please note the different energy scales of the subpanels. Positions 
ofq = (1/2 1/2 1/4) are over-plotted by dots. 

We note an almost linear increase of phonon energies with temperature. Interestingly, the 
dependence on Z point and q=(l/21/21/4)is very similar. For both g-points the calculated 
phonon energies are signihcantly too large. At q = (0 1/4 0) the temperature dependence 
is less pronounced and the agreement with the calculation is rather good. 

The effect of anharmonicity is investigated by frozen phonon calculations at q = (1/2 1/2 
1/4). The resulting phonon energy potential surfaces for different isotropic cell expansions 
and the corresponding displacement pattern are shown in Figure 6 together with quadratic 
(x^) and quartic (x^ and terms) hts. The potential energy curves show clear evidence 
of anharmonicity. All curves require x^ terms for a successful fit. A double-well potential 
forms for one percent expansion. 



















FIG. 4. TDS of Ti02 crystal, (a) and (b) Intensity distribution across high symmetry planes in 
reciprocal space, reconstructed from experiment (left part of individual panels), are compared to 
the calculation (right part of individual panels), (c) and (d) Iso-intensity distribution of TDS in 
3D around Q = (2.5 2.5 1) as obtained from experiment and calculation, respectively. The isovalue 
is chosen arbitrarily for best visualization. Colors indicate the distance from the (2.5 2.5 I) point. 
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FIG. 5. Temperature dependence of the phonon energies of the TA mode at selected q-values. (a) 
Energy transfers measured by IXS (connected blue symbols) and calculated values (black symbols 
at zero temperature) at q = (1/2 1/2 1/4), Z and q = (0 1/4 0). (b) Dispersion along T — Z, 
measured by IXS (connected symbols, 20K data at Z point only) and calculation (solid black line). 



O.BOr 


0.25 


0.20 

> 

0) 

0.15 

>> 

cn 


oi 

c 

LU 

0.10 


0.05 


0.00 


-0.05 


• • 0.0% isotropic expansion 
■ ■ 0 . 5 % 


- - quadratic fits 

— quartic fits 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 
Ti ion displacement (A) 



FIG. 6. Frozen-phonon potential-energy curves for different isotropic cell expansion (symbols®) for 
the soft TA mode at q = (1/2 1/2 1/4). The mode amplitude is given for Ti ion displacement. 
The dashed and full lines in the corresponding color show, respectively, quadratic and quartic fits 
to the DFT calculation. Eigenvectors are represented within the unit cell. 
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IV. DISCUSSION 


The combined use of the experimental techniques TDS and IXS enables to probe ex¬ 
tended regions in reciprocal space and thus allowed for a detailed measurement of the low 
energy TA phonon mode. A predicted dip structure around q = (1/2 1/2 1/4) could be 
conhrmed. It is worth noting, that this phenomenon was entirely missed by measuring the 
phonon dispersion along high symmetry directions. Modern techniques, including the cal¬ 
culations, allow scanning across the Brillouin zone and are thus sensitive to such anomalies. 
Experimental results on phonon energies and scattering intensities show that the harmonic 
DEPT calculation reproduces dispersion and eigenvectors of most phonon branches, includ¬ 
ing its rich variety of features. The shape of phonon energy surfaces is correctly described 
across the Brillouin zone in spite of significantly underestimating the phonon energy of the 
TA mode. 

The theoretical description of the lattice dynamics of Ti02 in the rutile structure from 
hrst principles methods is rather challenging. The mixed covalent bonding oi s — d orbitals 
is responsible for long-ranged force constants. In order to capture subtle features in the 
lattice dynamics, including the incipient ferroelectric behavior and the anomaly descried 
here, the calculations require careful convergence tests. The present calculation profits 
from accurate pseudopotentials and consequently precise wave functions providing lattice 
parameters within 0.5 - 1.0 % compared to the limit in LDA®. 

Some optical phonon branches, and most importantly to this study, the lowest energy 
TA mode are substantially underestimated in energy. Worthy of comment is the particular 
sensitivity of the TA phonon mode with anomalous low energy to convergence issues and 
external stress. An isotropic compression of the relaxed structure by 0.5 % brings the phonon 
frequencies closer to the low temperature experimental values in LDA calculations®. 

The temperature dependent experimental data and the frozen phonon calculations show 
the importance of anharmonic contributions to the lattice dynamics in Ti02. We hnd that 
the phonon energies at Z point and q = (1/2 1/2 1/4) increase almost linearly with tem¬ 
perature. Measurements of the phonon density of states as function of temperature by Lan 
et al.^ show that the entire TA branch stiffens with temperature. This observation is inline 
with our results. The change in energy at q = (0 1/4 0) is less pronounced and potentially 
provides important details how the phonon density of states broadens with temperature. 
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The mode level resolution of our measurements will allow for direct comparison with tem¬ 
perature dependent dispersion relations. A proper anharmonic treatment, which is beyond 
the scope of this work, is required for a correct description of the TA mode. Nevertheless, 
even without such treatment we can extract some information from the data presented in 
Figure 6. The phonon mode potential has strong positive fourth order components. This 
leads to an increase of the phonon frequency for large atomic displacements (i.e. at higher 
temperatures) which may explain at least part of the discrepancy between purely harmonic 
calculation in Figure 5b and non-zero temperature measurements in the same hgure. Fur¬ 
thermore, a rather strong dependence of fourth order components of the potential on unit 
cell volume indicates that thermal expansion is an important factor in the phonon frequency 
- temperature relationship of the investigated TA mode. 

A full quantum mechanical description of the lattice dynamics including anharmonic ef¬ 
fects would require the computation of the multi-dimensional potential-energy surface and 
solving Schrodingers equation for the relevant degrees of freedom. Ab initio molecular dy¬ 
namics simulation followed by a mapping onto a model Hamiltonian that describes the lattice 
dynamics^^ could be an option. This so called ’temperature dependent effective potential 
technique’ was successfully used to describe the temperature dependence of the phonon 
density of states in Ti 02 rutile® and phonon dispersions in VO 2 rutile^^, which is similarly 
affected by anharmonic contributions. The employed method should be reliable for classical 
nuclear motions but quantum effects need to be considered at low temperatures. Alterna¬ 
tive methods for studying anharmonicity have been proposed, for example in References'®'^®. 
The presented data provide important information for benchmarking such approaches. 

The harmonic DFPT calculation by Mitev et alA predicted the TA mode to become soft at 
q = (1/2 1/2 1/4) under isotropic expansion or uniaxial strain. Our temperature dependent 
measurements show the opposite under thermal expansion. This is caused by the fact that 
even at substantial (1%) expansion the energy difference between the global minimum and 
the local maximum in the potential is small in comparison with thermal energy. It suggests 
that the anharmonic contribution to the phonon energies at higher temperatures is much 
more important than the pure volume effect arising from thermal expansion. The open 
question to be addressed is what happens under application of strain. In fact, from the 
harmonic calculations we expect the TA mode to soften under application of external stress 
at small displacements. The application of uniaxial expansion could be realized by bending 
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thin films or using tensile load. TDS and IXS studies can be performed in such conditions 
in grazing incidence^® or transmission geometry, respectively. 

The experimental proof of the existence of the TA mode with an anomalous low energy 
minimum at q = (1/2 1/2 1/4) is important for the explanation of why perturbations on 
cleaved Ti02-surfaces have an effect deep into the bulk. Our results support the proposition 
of Mitev et al.^ for the anomalous slow and oscillatory convergence of thin him (’slab’) models 
in Ti02. The eigenvectors of the soft TA mode at q = (1/2 1/2 1/4) are very similar to layer 
displacements perpendicular to the (110) surface, see Figure 6. As the second layer below 
the surface has a structure similar to bulk Ti 02 ^°’^^, the high compliance of the eigenvectors 
allow for long-ranged propagation of surface displacements. The surface oxygen lone-pair 
orbitals are under-coordinated and result in strengthened Ti-0 bonding and shorter Ti-0 
distances between the topmost pair of planes and larger distances between this plane and 
the next lower plane. Depending whether a slab model consists of odd or even numbers of 
layers, the displacements are in or out of phase, respectively. 

From an experimental point of view this study is a benchmark for the sensitivity of diffuse 
x-ray scattering to particular phonon branches. In Ti02 the scattering from a specihc branch 
contributes up to 98% to the total TDS not only close to F. The shape of TDS can be rapidly 
explored over a large volume in reciprocal space with good momentum resolution. 


V. CONCLUSIONS 

The lattice dynamics of Ti02 exhibits a rich structure including an TA mode with anoma¬ 
lous low energy with a minimum at q = (1/2 1/2 1/4). With help of inelastic and diffuse 
x-ray scattering we were able to measure the energy of this mode over an extended region 
in reciprocal space and validate its eigenvectors by careful analysis of scattering intensi¬ 
ties. DFPT correctly described the shape of the dispersion relations including many subtle 
features. Absolute phonon energies are underestimated for several branches including the 
TA mode of interest. It shows a signihcant stress dependence and is thus very sensitive to 
the approximations used in the calculations. The temperature dependent study revealed 
a strong temperature variation of the soft TA mode over an extended region in reciprocal 
space which shows the importance of anharmonic contributions to the lattice dynamics in 
Ti02. Our experimental observation of the soft TA mode with a minimum at q = (1/2 1/2 
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1/4) finally supports the explication by Mitev et al.^ of anomalous convergence behavior in 
slab model calculations and real behavior in thin dims. 
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